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We discuss the stability of strangelets by considering dynamical chiral symmetry break- 
ing and confinement. We use a U(3)l X U(3)r symmetric Nambu-Jona-Lasinio model for 
chiral symmetry breaking supplemented by a boundary condition for confinement. It is 
' shown that strangelets with baryon number A<2xl0 3 can stably exist. For the observ- 

ables, we obtain the masses and the charge-to-baryon number ratios of the strangelets. 
These quantities are compared with the observed data of the exotic particles. 
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1 Introduction 



The strange matter, containing the u, d and s quarks, has been considered to be the ground 
state of QCD PJEIIHI, and expected to play an important role in the astrophysical phenomena 
such as the quark stars and the early universe. It is also interesting that droplets of the strange 
matter (strangelets) could be a candidate of the dark matter. There is also a possibility to 
observe strangelets in the relativistic heavy ion collisions. 

Whereas the physics of the quark matter is interesting in many respects, it is difficult to 
describe its properties directly from QCD. Until now, the stability of the strange matter has 
been discussed by using effective models of QCD. In the early stage, the MIT bag model was 
often used with an assumption that the strange matter could be treated as a system of free 
fermi gas in a bag 121 El E] There, the strange matter became stable than the ud quark 
matter due to the large number of the degrees of freedom by including the strangeness. In 
these works, not only the bulk quark matter, but also the strangelet of finite size has also been 
studied. It was then shown that the strangelets could be more stable than the normal nuclei. 

However, in the MIT bag model, an important feature of QCD is not incorporated, that is 
the dynamical chiral symmetry breaking. In fact, when it is taken into account, it has been 
shown that the strange matter cannot be absolutely stable |S1 13 El- The pattern of chiral 
symmetry restoration in the finite density quark matter is different for the ud quarks and the 
s quarks. The chiral symmetry for the ud quarks is sufficiently restored at stable densities 
{jib — 2 — 3n° B ), while that for the s quarks is still largely broken. Then, the transition from 
the ud quark matter to the strange matter by weak process is disfavored because of the large 
dynamical quark mass in the strange quark sector. This result is qualitatively very different 
from the results of the absolute stability obtained in the MIT bag model, in which only current 
quark masses are used. 

Though the strange matter of infinite volume was discussed by taking into account dynami- 
cal chiral symmetry breaking jSJIHllZj) the strangelet with finite volume is not discussed yet. In 
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this paper, we study the stability of the strangelet by considering dynamical chiral symmetry 
breaking. In section 2, we formulate the lagrangian to describe dynamical chiral symmetry 
breaking supplemented by the confinement which is treated approximately by a boundary con- 
dition. In section 3, numerical results are shown and we discuss the stability of the strangelet 
and present several observables. In section 4, we conclude our discussions. 



2 Formulation 

We consider the N JL model with a four-point quark interaction [H] . This interaction is respon- 
sible for the dynamical generation of the quark mass. For a finite system of strangelet, we 
also have to consider that quarks are confined in a cavity with finite volume. To incorporate 
the latter aspect, we introduce the boundary condition of the MIT bag model. Therefore, our 
model lagrangian for the strangelet is given by fOJ 



C 8 r 2 



^ a 11 



-tpM9(r-R)ij, 



where ip = (u,d,sY is the quark field and m° = diag(m°, m% m") current mass matrix. The 
second term of Eq. (JTJ) is the four-point quark interaction invariant under U(Nf) L x U(Nf) R 
symmetry, in which A a (a = 0, • • • , 8) are the Gell-Mann matrices normalized by trA a A 6 = 25 ab . 
In our formulation, we do not consider U(1)a breaking for simplicity. 

The last term in Eq. (0) has been used in the MIT bag model to impose quark confinement 
f2JE3E]- Assuming that the strangelet has a spherical shape, the step function 6(r — R) is 
introduced, where R is the bag radius. That term represents a quark mass term with M for 
the exterior region and quarks are confined in the region r<R by taking the limit M-^oo. It 
is well known that the last term of Eq. (pQ) breaks chiral symmetry explicitly at the bag surface. 
In order to recover chiral symmetry there, we need to introduce the chiral field (pion) which 
is coupled to the quarks at the bag boundary. This leads to the condition of the chiral bag 
model, in which the pion cloud exists around the bag, and the vacuum structure is modified in 
the bag due to the strong pion-quark coupling which is known as the chiral Casimir effect ^3] . 
Although, this property associated with the pion field does not appear explicitly in our model 
lagrangian Eq. (JIJ). the term of the NJL interaction in Eq. ([!} is considered to be responsible 
for the pion induced property. 

The parameters in Eq. (0), such as the coupling constant G, a three dimensional momentum 
cut-off A (see Eq. (JSJ below) and the current masses m" = m° d and m° s are determined to 
reproduce the pion mass m n = 0.139 GeV, the pion decay constant f n = 0.093 GeV and the 
averaged mass of the nucleon and the delta ttin+a = 1-134 GeV. The mass is considered to be 
related with the dynamical quark mass in the vacuum by mjv + A — 3m*. Then, we obtained the 
parameter set GA 2 = 4.7, A = 0.6 GeV and mJJ = m° = 5.9xl0~ 3 GeV. In this paper, we consider 
the current mass of the strange quark m° s as a free parameter. We show the results by setting 
m° = 0.1 GeV. Other choices of m° s do not affect our final conclusions qualitatively. 

Now, let us investigate chiral symmetry breaking in a quark bag. As usual, in the NJL 
interaction term in Eq. (fTj). we adopt the mean field approximation (qq) 2 — ► 2qq(qq) — (qq) 2 , 
where q = u,d and s, and solve the following gap equation 

m q = m° q -2G(qq). (2) 

Here we need to solve this equation for a quark bag. This requires a treatment of quark states in 
the quark bag with discretized energy levels, which is rather complicated. In order to simplify 
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the numerical calculations, we perform momentum integration with a density of states, which 
is approximately obtained by the multiple reflection expansion (MRE) jUE]- It is expressed by 
a smoothed function 



6tt 2 12tt 2 
Pmre{p, m, R) = 1 + -^fs(p/ m ) + 7^m2-fc(P/ TO ) + 



(3) 



where p is the momentum, m the dynamical quark mass and R the radius of the quark bag. 
The second and third terms in Eq. (J3J) are the correction terms by the surface and curvature 
effects. The functions fs and fc are given by 
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(4) 



In the limit m^O, and /c become constants; 

Ihn fs(jp/m) = 0, 

lim fc(p/m) = 

m^o JOW ; 24?r 2 



(5) 



By using the MRE method, the energy density e in a strangelet with a radius i? is given as 



E 

q=u,d,s 



4G 



' / „ Jp 2 + m 2 q p MRE (p,m q ,R) 



p 2 dp 



eo, 



(6) 



where the integral is modified by the density of state Eq. 0. Note that the dynamical quark 
mass m q in Eq. (JHJ is determined self-consistently by Eq. (j2J) as stated below. In Eq. (jOJ, 
f = N spin x A^" co ; or = 6 is the degrees of degeneracy of spin and color, and A in the integral 
is a three dimensional momentum cut-off. The value p q is the Fermi momentum which is 
determined by 



Pmre(p, m q ,R) 



p 2 dp 



n. 



(7) 



for a given quark number density n q for each flavor q = u,d and s. In Eq. (JHJ), the last term eo 
is the energy density in the chirally broken vacuum of infinite volume 



E 



q=u,d,s 



(m* — m° q ) 2 



AG 
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f 1 ^ I p 2 dp 

Jo V^ + < 2^ 



(8) 



where m* is the dynamical quark mass in the vacuum. Note that the energy density Eq. (jHJ) is 
written as a sum of the kinetic energy of the valence quarks and the effective bag constant just 
as in the MIT bag model, 



E - 

q =u,d,s 



p 2 + m\ p MRE {p, m q , R) C) 



p 2 dp 



effi 



where the effective bag constant is defined by 
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Note that the effective bag constant -B e // depends on the quark density, which is different from 
the constant B in the MIT bag model. Then, by taking de/dm q = 0, the gap equation Eq. (J2J) 
is written as; 



o f I p dp 

m <i = m q + 2Gu g^- J p\IP 2+m2 q PMRE(p,rn q ,R)-^-. (11) 



We mention that that the dynamical quark mass depends not only on the Fermi momentum, 
but also on the quark bag radius R. It is a characteristic feature of a finite size system. 

Now, by the energy density Eq. 0, the total energy of the strangelet with a radius R is 
given by 

E = eV + E c -^, (12) 

where V — (4n/3)R 3 is the volume of the strangelet. By assuming a uniform charge distribution 
in the strangelet, the Coulomb energy E c is given by 

where the total electric charge is given by Q = ^N u — ^N c i — ^N s with N q being the number of 
each quark q = u,d and s. The last term of Eq. (|12[) is a phenomenological zero point energy of 
the bag (a ~ 2.04). 

We obtain the energy of strangelet in the following way. First, we give a baryon number A 
and a strangeness fraction r s = N s /(N U + Nd + N s ) with N u = Nd. Then, for several radii R, we 
solve the gap equation Eq. (fTT|) and obtain the dynamical quark mass m q in the cavity, which 
is a function of the radius R. Then we find the minimum of the energy Eq. (|T2j) with respect 
to the radius R. 



3 Numerical result 

3.1 Chiral restoration in a cavity 

The confinement term in our model lagrangian Eq. (pQ) is responsible for the effects of the finite 
volume of the strangelet. In this subsection, we investigate dynamical chiral symmetry breaking 
in the empty cavity without valence quarks. 

In Fig. 21 we show the energy density Eq. (JHJ) as a function of the dynamical quark mass 
in the empty cavity with several radii for ud and s quarks (dashed lines). For comparison, 
the energy density of the bulk vacuum without the boundary condition is also plotted in the 
same figure (solid lines). Minimum points of the energy density provide the dynamical quark 
masses. Concerning the ud quark sector in the bulk vacuum, we obtain the dynamical quark 
mass m u = 0.378 GeV. On the other hand, in the cavity, the dynamical quark masses are 
m u = 0.322, 0.258 and 5.9 x 10- 3 (= m°) GeV for R = 20, 11.5 and 8 fm, respectively. We see 
that the dynamical quark mass becomes smaller as the radius decreases. This shows that chiral 
symmetry in the cavity tends to be restored. It is also true for the s quark sector. In the bulk 
vacuum, we obtain the dynamical quark mass m s = 0.539 GeV, while we find m s = 0.361, 0.175 
and 0.1(=m°) GeV for the radius R = 5.0, 3.1 and 2.5 fm, respectively. 

In Fig. the dynamical quark mass is shown as a function of the radius of the cavity. The 
chiral symmetry in the ud quark sector is restored at R— 11.5 fm, while that in the s quark 
sector is restored at _R = 3.1 fm. We see that the chiral restoration in the s quark is suppressed 



4 




m u [GeV] m s [GeV] 



Figure 1: The energy density e q as a function of the dynamical quark mass m q for various radius R. 
Left: u quark, right: s quark. 
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Figure 2: The dynamical quark mass m q of the u and s quarks as a function of the radius R. 



as compared with the ud quark. This is considered to be due to the large current mass of the s 
quark as compared with the ud quark. The difference in the tendency of the chiral symmetry 
restoration in the ud quark sector and the s quark sector causes the difference in the stability 
of the ud quark droplets and the strangelets. 

3.2 Stability of strangelet 

Now we discuss the stability of the strangelet. For this purpose, we add 3 A valence quarks in 
the cavity, and calculate the energy per baryon number E/A for several baryon numbers A and 
the strangeness fractions r s . In order to simplify the discussions, first, we fix the strangeness 
fraction r s = for the ud quark droplets and r s = 1/3 for the strangelets, respectively. In the 
case of r s = 1/3, the Coulomb energy of Eq. (JT3J) vanishes. In this subsection, we turn off the 
Coulomb term. 

The cavity radius R of a strangelet is determined by the variation of E/A. In Fig. Efa), we 
show the energy per baryon number E/A as a function of the cavity radius R for the baryon 
numbers A = 10 2 , 10 3 and 10 4 . The minimum of E/A gives the energy and the radius of the 
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Figure 3: Left: (a) The energy per baryon number E/A as a function of the radius R for the ud quark 
droplet (r s = 0) and the strangelet (r s = l/3). Right: (b) The dynamical quark mass m q of the ud and 
s quarks in the strangelets with r s = l/3. 

strangelet. The existence of the minimum is understood in the following way. In the total 
energy Eq. ()12j) and (|5)i. there are two terms: the kinetic energy from the valence quarks and 
the volume energy B e ffV from the effective bag constant. As the radius becomes large, the 
kinetic energy decreases, while the volume energy increases. On the other hand, as the radius 
becomes small, the kinetic energy increases, while the volume energy decreases. Then, we find 
the equilibrium radius balanced by the kinetic energy and the volume energy. 

For example, for the baryon number A = 10 4 , the energy per baryon number in the ud quark 
droplet and the strangelet are (E/A) ud — 1.27 GeV and (E/A) uds = 1.33 GeV, respectively. For 
the baryon number A = 10 2 , the energy per baryon number in the ud quark droplet and the 
strangelet are {E/A) ud = 1.60 GeV and (E/A) uds = 1.48 GeV, respectively. In our results, the 
strangelets are more stable than the ud quark droplets for smaller baryon numbers A<2xl0 3 . 
The stability of the strangelets with small baryon numbers is very much different from the result 
for the bulk quark matter, where the strange matter with infinite volume is not absolutely stable 
O 13 E] • This is because of the effect of the confinement leading to the restoration of chiral 
symmetry in the cavity. In order to show the restoration of chiral symmetry in the strangelets 
in Fig. 0(b), the dynamical quark masses m u and m s of the ud and s quarks are shown as 
functions of the cavity radius R. We see that chiral symmetry of the ud quark and s quark in 
the strangelet has a tendency to be restored for small radii as seen in the empty cavity. 

In Fig. IHa), we show explicitly the energy per baryon number E/A as a function of the 
baryon number A for ud quark droplets and strangelets. It is shown that the strangelets are 
more stable than the ud quark droplets for the baryon number A < 2 x 10 3 . It is generally 
expected that the strange matter can be more stable than the ud quark matter, when the 
dynamical quark mass of s quark m s is smaller than the Fermi energy ep )U of the ud quark. 
When this relation is satisfied, the weak transition from ud quarks to s quarks can occur by 
the weak processes u— >d+e + +v e and u+d— >u+s. In Fig. 0(b), we compare m s and ep >u in the 
ud quark droplets. It is shown that m s <e j p jU is satisfied in the quark droplets with v4<2xl0 3 . 

So far, we have discussed the stability of the strangelets as compared with the ud quark 
droplets. We have found that the absolute values of the energy of the strangelet for v4<2xl0 3 
is larger than the masses of normal nuclei. Consequently, the strangelets may not be absolutely 
stable as the QCD matter. However, if there are many s quarks in the strangelets, the decay of 
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Figure 4: Left: (a) The energy per baryon number E/A as a function of the baryon number A. The 
solid line for the strangelets (r s = l/3) and the dashed line for the ud quark droplets (r s = 0). Right: 
(b) The dynamical quark mass m u and m s of the u and s quarks, and the Fermi energy €f, u of the u 
quark in the ud quark droplets with r s = 0. 

the strangelets to the normal nuclei would take a long time due to the many steps in the weak 
processes. Therefore the strangelets could survive as stable particles, once they are formed. 

3.3 Observables 

When we consider that the strangelets are formed in the QCD phase transition in the early 
universe and/or in the explosions of the strange stars, the remaining strangelets could be 
observed in the cosmic rays. They would be observed as exotic particles with a large mass and 
a small electric charge. In order to identify such heavy particles, the charge-to-baryon number 
ratio is an important quantity. For this purpose, let us concentrate on the stable strangelets of 
A<2xl0 3 . 

In the previous subsection, we have fixed the ratio of s quarks to be r s = and 1/3. In 
the following, we consider a variation of the energy Eq. (|T^j) with respect to r s . The resulting 
strangeness fraction r s is plotted as a function of the baryon number A in Fig. 03 We see that 
the strangeness fraction of the strangelets with A<2x 10 3 is close to r s = l/3. 

Once the strangeness fraction r s is obtained, the number of quarks N q for q = u,d and 
s are also obtained. Then, the electric charge Q of the strangelets are calculated from Q = 
^N u — ^Nd—^N s . In Fig. El we show also the charge-to-baryon number ratio Q/A as a function 
of the baryon number A. The electric charge of the strangelets are of order a few percents of the 
normal nuclei. Such strangelets can have large baryon number, since the Coulomb instability 
can be negligible. Thus, our results show that the strangelets would be the exotic particles 
with small charge-to-baryon number ratio as compared with the normal nuclei. 

Let us compare our theoretical results and the existing data which was reported in the 
observation of exotic particles in the cosmic rays. We show the baryon number and the charge- 
to-baryon number ratio of the observed particles in Table ^ The baryon numbers A of these 
exotic particles are from the order of hundreds to thousands. The charge-to-baryon number 
Q/A is around 0.04. It is interesting that these observed values are very close to our theoretical 
results. These exotic particles could be candidate of strangelets. 

In Fig. El we show the radius R of the strangelets as a function of the baryon number A. 
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Figure 5: The strangeness fraction r s and the Figure 6: The radius R of the strangelet as a 
electric charge Q/A of the strangelets as a func- function of the baryon number A. 
tion of the baryon number A. 

The relation between R and A is expressed approximately by 

R = r A 1/3 , (14) 

with r — 0.57 fm. Then, the baryon number density in the strangelet is estimated as % = 
A/ (AnR 3 /3) = 7.6 n B with the normal nuclear matter density n° B = 0.17 fm -3 . Recently, the 
K-nuclei including the strangeness is discussed to have such a high density (2H1 EE] . Though 
the relation between the strangelets and the X-nuclei is not yet clear, it is interesting that such 
high baryon number density in the strangelets is comparable with the result in the ^-nuclei. 

So far, we did not include electrons in our discussions. We show that the electrons play only 
minor role in the strangelets. From Fig. the electric charge of the strangelets are at most 
Q~55 for A = 2x 10 3 , due to small Q/t4~0.026. Assuming an electron around the strangelets 
with such electric charge, the de Broglie wave length of the electron is estimated by the energy 
of the electron E e = ^Jp 2 + m 2 — e 2 Q/r, where p is the momentum of the electron, m e = 0.51 
MeV the mass of the electron and r the distance from the strangelet. Then, we find that the de 
Broglie wave length of the electron is A~960 fm for A = 2x 10 3 . That is much longer than the 
radius of the strangelet, which is shown in Fig. El Therefore, the electrons exist almost outside 
of the strangelet. If the size of the strangelets could be larger than the electron de Broglie wave 
length, we need to consider the electrons in the strangelets. 



baryon number A 


charge-to-baryon number Q/A 




A ~ 350 -450 


0.03-0.04 


M 


A~460 


0.043 


!2H 


A > 1000 


0.046 


[23 


A -370 


0.038 


[23,121 



Table 1: The baryon number A and the charge-to-baryon number Q/A from cosmic ray experiments. 
See also [T§] . 
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4 Conclusion 



We discussed the stability of the strangelet by considering dynamical chiral symmetry breaking. 
We used the NJL interaction for chiral symmetry breaking. In addition, we incorporated the 
spherical cavity for the confinement of quarks by the boundary condition approximately. In 
the mean field approximation in the finite volume system, we obtained the gap equation for 
the dynamical generation of the quark mass. Then, we obtained the energy of the strangelets. 

As a result, it was shown that chiral symmetry tended to be restored in the cavity at small 
radii. Then, the dynamical quark mass became small as compared with that in the vacuum of 
infinite volume, and the strange quark mass can be smaller than the Fermi energy of the ud 
quarks in the droplets. We investigated the stability of the strangelet for several baryon number 
A and the strangeness fraction r s = and 1/3. Then, it was shown that the strangelets are more 
stable than the ud quark droplets for the baryon number A < 2 x 10 3 for the s quark current 
mass m° = 0.1 GeV. Our result did not change qualitatively for the case of m° = 0.18 GeV, in 
which we obtained the stable strangelets for A<0.5x 10 3 . We obtained the charge-to-baryon 
number of the strangelets, which is consistent with the data which was reported to be observed 
in the heavy particles in the cosmic rays. 

Finally, we comment on the further developments in the study of the strangelets. The 
energy per baryon number E/A in the stable strangelets with A<2x 10 3 seems to have larger 
energy than ordinary nuclei. For the more realistic discussions, the color exchange interaction 
between dynamical quarks should be introduced. It is known in the quark model that the 
color exchange interaction splits the mass of the nucleon N and the delta A, while, in our 
calculation, we considered that the masses of the N and A are degenerate. Furthermore, the 
quark-quark pairing, namely the color superconductivity, in the quark droplets may affect the 
energy of the strangelet [TBI HZl HHj- Especially, it is discussed that the strange matter can 
form the color-flavor locked phase. Such a condensed state could play an important role in the 
strangelets. 

We would like to express our gratitude to Prof. S. Raha and Prof. N. Sandulescu for inter- 
esting suggestions to start our study of strangelets. We would express thanks to Prof. Itahashi 
for fruitful suggestions from the view of experiments. We wish to thank Dr. O. Kiriyama for 
many discussions from an early stage of our study. 
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